Ef»A/60Q/A-9?/058 





PB97-195499 


RATION OF SATELLITE DATA IN REGIONAL AIR QUALITY 


Richard T. McNider, William B. Norris, and Daniel Casey 

Earth System Science Laboratory 
University of Alabama in Huntsville 
Huntsville, Alabama 


Jonathan E. Pleim and Shawn J. Roselle 

Atmospheric Science Modeling Division 

Air Resources Laboratory 

National Oceanic and Atmospheric Administration 

Research Triangle Park, North Carolina 

(on assignment to the National Exposure Research Laboratory, U.S. EPA) 


William M. Lapenta 

NASA Marshall Space Flight Center 
Global Hydrology and Climate Center 
Huntsville, Alabama 


INTRODUCTION 


In terms of important uncertainty in regional-scale air-pollution models, probably no 
other aspect ranks any higher than the current ability to specify clouds and soil moisture on 
the regional scale. Because clouds in models are highly parameterized, the ability of models 
to predict the correct spatial and radiative characteristics is highly suspect and subject to large 
error. While considerable advances have been made in the assimilation of winds and 
temperatures into regional models (Stauffer and Seaman, 1990), the poor representation of 
cloud fields from point measurements at National Weather Service stations and the almost 
total absence of surface moisture availability observations has made assimilation of these 
varmbles difficult to impossible. Yet, the correct inclusion of clouds and surface moisture are 
°t rirst-order importance in regional-scale photochemistry. Consider the following Points 
relative to these variables. & y 
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OTKaoris models has gone beyond our ability to specify surface temperatures m a sbrtS 

LgM 

or aciwdy transpiring vegetation can sharply reduce temperatures over that of di? surfaces 

wKX i,y ^ ^’to^Wcrrorr^ 

back into the emission errors discussed above. reeaing 

he jg hts .. in regional-scale models, which affect pollution and orecursor 
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to summary, the above discussion shows the importance of cloud* and surface mnictntv* 
availability , n regional models. While considerable emphasis hTbin plSoTS 
f’ r f!! on ® regional scale models, the winds affect the distribution of pollutants but clouds 
!*. ™ 0,st V re a f ^ aiIabl ! It y highly affect the photochemical production. (The clouds and 
moisture also > affect wmds-see below). The domino effect of clouds and rnoistum throSh 

models. S,S ra S ’ emisslons ’ m,xm 8 heights, etc., makes them a pivotal element in regional 

„ . ^ is the purpose of this paper to describe methods of satellite remote sensing that can be 

oh l to Sp - Clfy > C 0l I ds o nd sarface rnoisture in photochemical models with improved fidelity 
on he regional scale Specifically, geostationary data is used because of STmtfoSl aSd 


INSOLATION 


.. im T j lc net solar radiation at the surface provides the prime source of energy controlling the 
umal variation in temperature in the surface energy budget in regional S scale models. In 



f 5 tore > ^ " et Nation at the surface is determined by the surface 
albedo, reflection, and absorption of solar radiation by the the clear atmn<nhp~> t h* 
reflection and absorption by the cloudy atmosphere (including aerosols). ^ 

Surface Albedo 

aIbe< ?, 0 “ regional-scale models such as MM5 (Grell et al , 1994) or RAMS 
1S u ^ uall y estimated using land-use types. However, the relationship between 
gross land-use characteristics and its radiative properties is not always well defined 

5f!lJ? 0n f Uy albed ? c f n ° ba ”S e due to meteorological conditions tod anthropogenic 

W ^- h T T ? ncluded to T6v ^ AQ ^-nst data bases. Satellites 
on the other hand, provide a direct radiative measurement of reflected radiation although 

1980) rctatl ° n dUC t0 Vl6W ang e ^ bi ' directional effects must be considered (Gautier et al., 

The geostationary series of satellites (GOES) operated by NOAA in the li s over the 
« wo decades provide satellite coverage over the ifs. thatSt. te^ed it, mSet 
a n ® w vers . 10 { 1 of geostationary satellites beginning with GOES-8 was launched with 

Ss n -Ri C «S-f nStl ? ** P revious G0ES ( ie ? G0ES ? and earlier). The GOES-7 
( OES 8) satellite returns the magnitude of upwcllitig radiance in the visible as brightness 
coante to the rsrige 0-63 (0-1023). The lowest Lnts Sise wheVStdTor no ?Cd U fS 
and the reflation is primarily, if not entirely, from the earth’s surface. Using the physical 
retrieval method of Gautier et al. (1980) and Diak and Gautier (1983), the surface albedos 
can be recovered from the brightness counts returned by GOES satellites. EecS 
requires hourly surface albedos obtained from clear-sky brightness counts. If a single cloud- 

dS‘f^ aVa, » blef ° re! 1' hour '*,*&&’ *« brightness counts could to obS 
directly from them. However, because cloudy skies are so common, a single, cloud-free 
image is usually not available. Experience has shown that for a given daylight hour images 
h^f.r a ? 2 ?' 3 .9 da £ s ar ® need f d . t0 obtain a stable minimum brightness count for that 

• thB S ° U ? e f s - t dmm i Summer when cumulus clouds are ubiquitous. 
Brightness counts are converted into reflectances using a calibration curve uniaue to each 
sateUite. This approach has the inherent ability to account for both spatial and temporal 

— d ^ e t0 s6ll . t .yp e ’ vegetation, tod time of day and year. Figure 1 shows a 
clear-sky albedo derived using this technique tod converted to to 80-km grid used in MM5. 



Figure 1. Clear-sky albedo using minimum brightness from GOES-7 
an 80-km MM5 domain. 


images du-ing July 1988 gridded to 


Cloud Albedo 


• ° nc ®. toe dejtf-sky albedo and the brightness count are known for a given satellite-image 

I)ixel surface insolation at the pixel level (usually 1, 4, or 8 km) can be calculS from &e 

mqsfTh radia S? :tranSfer m<?de ! 0f , Ga 1 U , tier et al (I980) 38 Ascribed in McNider etal. 

TTu S a single cloud H yef ‘ Abovfe cloud J ayer, radiation is Rayleigh 
!2S?i abs ff by water vapor. In the cloud layer, radiation is scattered a5d 
abs 9£ b ® d, ® e ^ ow cloud layer, radiation is absorbed by water vapor. For the scattering 
coeflreients, we use the parameterization originally presented by Kondratyev and 

modified by j Atwater and Brown (1974). For the water-vapor absorption Coefficients, we use 
an empincal formulation of MacDonald (1960). For in-cloud absoiption we use a step 
function that depends on brightness count (McNider et al, 1995). This results in a quadratic 
equation m cloud albedo. Once known the cloud albedo can be used to calculate downwelline 
solar radiation and insolation at the surface. * 

A flow chart Of the computational procedure is shown in Figure 2. The procedure yields 
surface insolation at each pixel in an image. Such images can be gridded and values tor all 
pixels within a gnd cell averaged to produce hourly input fields for assimilation into 
photochemical models or meteorological models. Figure 3 shows the surface insolations 
computed by this technique for Julian day 216 on an 80-km MM5 grid covering most of the 



Figure 2. Flow chart of computational procedures for determining surface insolation. 
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Figure 3. Surface insolation computed from a Satellite image for an 80-km grid for 1800 GMT Julian Day 


While the GOES visible sensors do not have on-board calibration, prelaunch calibration 
curves (Raphael and Hay, 1984) and ground-truth data (Taipley, 1979) can approximate the 
needed calibration information. To provide some understanding of the fidelity of die 
technique a comparison of observed versus satellite derived insolation values for the period 
July 31-August 7, 1988 was made. Observed data were taken from available surface 
pyranometer archives (National Renewable Energy Laboratory, 1992). Because the GOES 
sensor measures a near instantaneous value of clouds and the pyranometer data are hourly 
averages exact agreement cannot be expected. Additionally, the absolute navigation of the 
GOES-7 image is of order 10 km. Thus, in making the comparison between the satellite data 
and pyranometer data the best-fit satellite pixel value within 20 km of the pyranometer site 
was utilized. Figure 4 shows this comparison. While Scatter exists, the key fact is that die 
calibration data used shows almost no bias in the insolation values. Because the insolation 
calculations are based on rather straightforward radiative principals and insolation is fairly 
insensitive to uncertain specifications such as cloud absorption (McNider etal., 1995), it is 
felt that absolute and especially the relative spatial accuracy is quite high. 


Hourly Surface Insolation 



Figure 4. Comparison between satellite-derived and surface pyranometer-measured insolation values for the 
period July 31-August 8, 1988. The best-fit pixel was used within a 20-km radius of the observation site. 
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PHOTOLYSIS RATES 

Photochemical modeling systems differ in their sources of photodissociation constants 
In some cases these are obtained from radiative transfer models that include the scattering 
?°sorbm g , or reflecting properties of atmospheric aerosols and clouds (<?.*., Ruggaber et al 
1994; Rosselle etal, 1995). In other cases radiative transfer models assume clear skies * and 
correction factors are later applied to the results to obtain the desired cloudy-sky values (e.g. 
Chang et al, 1987). In either approach cloud information is essential. ’ 

^ er , t0 ^ mput * clou< * effects on photolysis rates three fundamental parameters are 
needed: (1) Cloud-layer transmittance, (2) cloud-top elevation, and (3) cloud-layer bottom 
elevation. In tire past in photochemical models these parameters Were estimated based on 
hydrological information from the meteorological processor. Often the information used to 
recover these parameters is fairly indirect. For example, RADM approximates the cloud 
optical depth f with the parameterization 


T = iL con^ctd/ 2 PH 2 0 r 

where is the mean condensed water content, Az cW is the mean depth of the cloud layer, 
r is the mean cloud drop radius, and p/^is the density of water (Chang et al ., 1987). In 
RADM, constant values are assumed for each factor except &z dd , which is obtained from the 
meteorological model. This cloud optical depth is then used in a Beer’s Law type of 
formulation to estimate cloud transmittance. In practice the cloud top is often determined from 
relative humidity thresholds and cloud base from computed lifting condensation levels We 
now propose techniques that use satellite data to estimate two of the needed parameters- 
cloud transmittance and cloud-top elevations and discuss options for determining cloud base. 



Figure 5. Satellite-derived transmittance field for 1800 GMT Julian Day 216 1988. 


Satellite-Derived Transmittance 

Cloud transmittance is fundamentally a radiative property and can be recovered from the 
requirement that cloud transmittance, reflection and absorption sum to unity. We make the 
assumption that the broad-band transmittance determined from the satellite sensor (0;52 - 
0.72 pm) is the approximate transmittance required in the ultimate photolysis calculations. 
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Since cloud reflectance 1 $ determined above in the insolation step and if the same absorptivity 
function is used then cloud transmittance can be computed. Absorptivity is the least well 
known parameter in the present scheme. McNider et al (1995) showed in sensitivity tests 
matevert making extreme changes in the absorptivity function Changes the tranmittance only 
slightly. Figure 5 shows the transmittance field determined by Satellite for a RADM 80-km 
grid using these techniques described above. 

Satellite-Derived Cloud-fop Elevations 

Cloud-top elevations Can be estimated using infrared satellite images from atmospheric 
clear-window channels. Just as with visible images, GOES-8 returns infrared measurements 
m discrete counts in the range 0-1023. These values correspond to blackbody temperatures 
and m clear-window channels represent the temperature of the Earth's Surface in clear regions 
and cloud-top temperatures in cloddy regions. If a temperature sounding is available near the 
location of a pixel at the time of the satellite measurement, the pressure level, and hence the 
elevation, at cloud top can be determined. At this point in the development of the technique 
for a given pixel we are using the grid point sounding from the meteorological model. Figure 
6 shows satellite-derived cloud top elevations. 



Figure 6 . Satellite-derived cloud-top elevations for 1800 GMT Julian Day 216 1988. 


Cloud-Base Elevations 

GOES satellites cannot provide any guidance on the cloud base since infrared techniques 
determine the depth of clouds. However, auxiliary thermodynamic structure 
information can be used to compute the lifting condensation level (LCL) which gives a lower 
limit on the cloud-base elevation. This auxiliary information can come from prognostic 
mesoscale models or rawinsonde observations. 

Satellite-Derived Photolysis Fields 

The cloud characteristics defined above were utilized in a photolysis model inec porated 
in a new version of the RADM regional photochemical model. Figures 7 and 8 show a 
comparison between the photolysis fields generated from the satellite-derived cloud 
characteristics and the diagnostic determination of cloud characteristics from the prognostic 
meteorological model fields. The differences are quite dramatic for this hour, with the 



diagnostic models overestimating cloud coverage. Analysis of this case is continuing but it 
perhaps points to a major problem in that the highly parameterized cloud characteristics mav 
not be capturing the actual cloud fields with fidelity. y 



Figure 7. Satellite-derived photolysis fields (J N02 ) 1800 GMT Julian Day 216, 1988. 



Figure 8. Photolysis fields (J N02 ) derived from prognostic model output 1800 GMT Julian Day 216 1988. 


SURFACE MOISTURE AVAILABILITY 

As mentioned in the introduction, biogenic and soil NO, emissions rates highly depend 
on surface temperature. Mixing heights and turbulent dispersion characteristics are also 
dependent upon surface temperatures. Surface temperatures are highly sensitive to 
specification of surface moisture (Pleim and Xiu, 1995), yet moisture is not readily available 
as an observable and must be estimated using vegetation and soil moisture parametefizations 
that contain difficult-to-specify parameters such as deep root-zone moisture, stomata! 
conductance. To avoid these problems we have developed a method employing infrared 
satellite images that is capable of assimilating satellite-observed surface temperature rates of 
change into boundary-layer models in a thermodynamically consistent manner that recovers 
soil moisture. The technique requires the use of GOES-8 derived infrared images, land 
surface temperatures (LSTs), surface albedo, arid insolation over the time period of interest 
It is based upon adjusting the model’s bulk moisture availability so that the model’s rate of 
change of LST agrees more closely with the “observed” rate of change of LST as derived 
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from the satellite. A critical assumption is that the availability of moisture (either from the 
soil or vegetation) during mid-moming hours is die least known term in the model’s surface 
SSy . bud *« (Weoel ft <21 1984; Carlson 1986). Thoreforo, fte slmutattd EhMfla? 
which is a function of surface moisture availability, is adjusted based upon differences 
between the modeled and satellite-derived LST tendencies. 

r , w ® P resen f/ esuIt ? fr° m the assimilation technique as applied within MM5 for 7 
mly 1995 using a grid resolution of 12 km. The region of interest was over Oklahoma and 
Kansas where a west-east vegetation gradient existed (wheat stubble to the west, deciduous 
forest to the east). Given this variation in vegetation, one would expect the satellite-observed 
surface thermal response to differ across the region. The LST is derived via a physical split 

r, nd ?^^?^ Guilloi y etal ' 1993 ) from GOES 8 data from 1200 to 2300 UTC on 7 
July 1995. The LST tendencies are fairly uniform across the region from 1200 to 1400 UTC 
However, between 1400 and 1700 UTC the LST tendency in eastern Oklahoma (5°C per 3 h) 
is considerably lower than that farther West (10 4 C per 3 h). This spatial variation is consistent 
with the west-east vegetation gradient that existed across the region. 

TCie assimilation scheme in MM5 adjusts the evaporative flux by altering the moisture 
availability parameter, M, which represents the fraction of possible evaporation for a 
saturated surface. When the Standard procedure of specifying M as a function of land-use 
category was used to initialize the model, the field exhibits an unrealistic, blocky structure 
(Figure 9, left panel). J 



Figure 9. Model moisture availability, M, representing the fraction of possible evaporation for a saturated 
surface. The left panel (a) shows the moisture availability based upon land-use categories while the right 
panel (b) displays it after assimilating GOES data. 


After assimilating GOES LST tendencies from 1400 UTC to 1800 UTC (900 to 1300 
local standard time), the adjusted M exhibits a more realistic structure (Figure 9, right panel) 
which matched well with Normalized Difference Vegetation Index (NDVI) (not shown). A 
more stringent check of the technique was done by comparing modeled surface evaporative 
flux data from observations obtained via Energy Balance Bowen Ratio (EBBR) systems 
deployed at the Southern Great Plains (SGP) Cloud and Radiation Testbed (CART) sites for 
7jfnH «,-? epart . m6nt ? f . Energy’s Atmospheric Radiation (ARM) Program (Splitt et ai, 
tiV ^*2 ° Ut t ie assim ilation, MM5 greatly overestimates the daytime evaporative flux by 
300 Wm (150%) on both 7 and 8 July. However, use of the satellite data produces a much 
more realistic time series with an average difference for the 36h period of less than 40 Wm' 2 . 





10 


REFERENCES 

Atwater, M. A. and P. S. Brown, Jr., 1974, Numerical calculation of the latitudinal variation of solar 
radiation for an atmosphere of Varying opacity, J. Appl. Meteor., 13:289-297. 

Carlson, T. N 1986, Regional scale estimates of surface moisture availability and thermal inertia using 
remote thermal measurements. Remote Sensing Rev., 1: 197-246. 

Carlson, T. N., J. K. Dodd, S. G. Benjamin, and J. N. Cooper, 1981, Satellite estimation of the surface 
energy balance, moisture availability and thermal inertia, J. Appl. Meteor., 20:67-87. 

Chang, J. S..R. A. Brost, L S. A. Isaksen, S. Madronich, P. Middleton, W. R. Stockwell, and C. J. 

Walcek, 1987, A three-dimensional eulerian acid deposition model: physical Concepts and formulation 
J. Geophys. ReS., 92:D12, 14681-14700. 

Deardorff, J.,1974, Three-dimensional numerical study of the height and mean structure of the planetary 
boundary layer, Bound-layer Meteor. 15:1241-1251. 

Diak, G. R. and C. Gautier, 1983, Improvements to a simple physical model for estimating insolation from 
GOBS data, J. Appl. M"eor., 22:505-508. 

Dunker, A. M., 1980, The response of an atmospheric re/action-transport model to changes in input function 
Amos. Environ., 14:671-679. 

Gautier, C., G. Diak, and S. Masse, 1980, A simple physical model to estimate incident solar radiation at the 
surface from GOES satellite data, J. Appl. Meteor., 19:1005-1012. 

Grell, G. A., J. Dudhia, and D. R. Stauffer, 1994, ,4 Description of the Fifth-Generation Penn State/NCAR 
Mesoscale Model (MM5). NCAR Technical Note NCAR/TN-398+STR, Nationial Center for 
Atmospheric Research, Boulder, Colorado. 

Guillory, A. R., G. J. Jcdlovec, and H. E. Fuelberg, 1993, A technique for deriving column-ihtegrated water 
content using VAS split-window data, J. Appl. Meteor,, 32:1226-1241. 

Kondratyev, K. Y., 1969, Radiation in the Atmosphere., Academic Press, New York. 

McDonald, J. E., 1960, Direct absorption of solar radiation by atmospheric water Vapor, J. of Meteor., 
17:319-328. 

McNider, R. T., J. A. Song, and S. Q. Kidder, 1995, Assimilation Of GOES-derived solar insolation into a 
mesoscale model for studies of cloud shading effects, Int. J. Remote Sensing, 16:2207-2231 

National Renewable Energy Laboratory, 1992, National Solar Radiation Data Base User's Manual (1961- 
1990), Golden, Colorado. 

Pielke R. A., W. R. Cotton, R. L. Walko, C J. Tremback, W. A. Lyons, L. D. Grasso, M. E. Nicholls, M. 
D, Moran, D. A. Wesley, T. J. Lee, and J. H. Copeland, 1992, A compreheriseive meteorological 
modeling $ysteiii~RAMS, Meteor. Atmos . Phys 49:69-91. 

Pleim, J. E. and A. Xiu, 1995, Development and testing of a surface flux and planetary boundary layer model 
for application in mesoscale models. /. Appl Meteor 34:16-32. 

Raphael, C., and Hay, J. E., 1984, Ah assessment of models which use satellite dta to estimate solar 
irradiance at the earth’s surface, J. Climate and Appl. Meteor „ 23:832-844. 

Roselle, S. L, A. F. Hanna, Y. Lu, J. C. Jang, K. L. Scherc, J, E. Pleim, 1995, Refined photolysis rates for 
advanced air quality modeling systems, in Proceedings of the A&WMA Conference on the Applications 
of Air Pollution Meteorology. 

Ruggabef, A., R. Dlugi, and T. Nakajima, 1994, Modelling radiation quantities and photolysis frequencies in 
the troposphere, J. Atmos. Chetn ., , 18:171-210. 

Seinfeld, J. H., 1988, Ozone Air Quality Models: A critical review, J. Air Poll. Control Assoc,. 38:616- 
645. 

Splitt, M. E. and D. L. Sisterson, 1995, Site Scientific Mission Plan for the Southern Great Plains CART 
Site: July-December 1995, ARM-95-002, Argonne National Laboratory, Argonne, Illinois. 

Stauffer, D. R, and N. L. Seaman, 1990, Use of four-dimensional data assimilation in a limited-area 

mesoscale model. Part I: Experiments with synoptic-scale data, Mon. Wea. Rev., liS: 1250-1277. 

Tarpley, J. D M 1979, Estimating incident solar radiation at the surface from geostationary satellite data J of 
Appl, Meteor., 18:1172-1181. * J 

Tingey, D. T., M. Manning, L. C. Grothaus, and W. F« Burns, 1979, The influence of light and temperature 
on isoprene emission rates from live oak, Physiol. Plant , 47:112-118. 

Wetzel, P. J., 1984, Determining soil moisture from geosynchronous satellite infrared data* A feasibility 
study, J. Climate Appl. Meteor ., 23:375-391. 

Zimmerman, P. R., J. P. Greenberg, and C. E. Wcstberg, 1988, Measurements of atmospheric hydrocarbons 
and biogenic emission fluxes in the Amazon boundary layer, J-. Geophys. Res., 93:1 407-1 41 6. 




TECHNICAL REPORT DATA 


1, REPORT' NO. 2. 

EPA/600/A-37/Q58 


4. TittiE AND SUBTITLE 

Assimilation of Satellite Sata in Regional Air 
Quality Models 


5. REPORT DATE 


6. PERFORMING ORGANIZATION CODE 


7. AUTHOR (S) 

R. T* McNider 1 , W.B . Nbrris 1 , D. Casey 1 , J.E. Pleim 2 , 

S. J. RoSelle 2 , and K,M. Lapenta 3 


9. PERFORMING organisation name and address 

^arth System Science Laboratory, University of 
Alabama-Huntsville, Huntsville, AL 
2 Same as Block 12 

3 NASA Marshall Space Flight Center, Global Hydrology 
& Climate Center, Huntsville, AL 


12. SPONSORING AGENCY NAME AND ADDRESS 

U.S . Environmental Protection Agency 
Office of Research and Development 
National Exposure Research Laboratory 
Research Triangle Park, NC 27711 


8. PERFORMING ORGANIZATION REPORT NO. 


10. PROGRAM ELEMENT NO. 


11. CONTRACT/GRANT NO. 


13. TYPE OF REPORT AND PERIOD COVERED 


Proceedings, FY-97 


14. SPONSORING AGENCY CODE 


EPA/ 600/9 



16. ABSTRACT 

In terms of important uncertainty in regional -scale air-pollution models, probably 
no other aspect ranks any higher than the current ability to specify clouds and soil 
moisture on the regional scale. Because clouds in models are highly parameterized, 
the ability of models to predict the correct spatial and radiative characteristics 
is highly suspect and subject to large error. While considerable advances have been 
made in the assimilation of winds and temperatures into regional models (Stauffer 
and Seaman, 1990) , the poor representation of cloud fields from point measurements 
at National Weather Service stations and the almost total absence of surface 
moisture availability observations has made assimilation of these variables 
difficult to impossible. Yet, the correct inclusion of clouds and surface moisture 
are of first-order importance in regional-scale photochemistry. 


KEY WORDS AND DOCUMENT ANALYSIS 


b. IDENTIFIERS/ OPEN ENDED 
TERMS . . * 



18. DISTRIBUTION STATEMENT 



19. SECURITY CLASS (This 
Report) 

UNCLASSIFIED 

29. SECURITY CLASS (This 
Page; 

UNCLASSIFIED 


21. NO. OF PAGES 


22. PRICE 


r.Mr €^p\ \C©OAq. s S'.rvv,! 
















« 0) 0 
O£o£ 

& v% v 

gli-2 

SIee 

$ o S® 

«4- a 
C id e-»+s 

U -■ <4- y_ 
m M ,w ^ 

4-* p wTu 

8 |«o 

S § 0 c 
is® o 

C«3<3 
ouo; 
Q.JS >c 

4- d 0) 0 
0 ££§ 
. g T? 
OjC® 

$!$? 

2o.2.£ 


Reproduced by NTIS 

National Technical Information Service 
Springfieid, VA 22161 


This report was printed specifically for your order 
from nearfy 3 million titles available in our collection. 

For economy and efficiency, NTIS does not maintain stock of its vast 
collection of technical reports. Rather, most documents are printed for 
each order. Documents that are not in electronic format are reproduced 
from master archival copies and aft the best possible reproductions 
available. If you have any questions concerning tills document or any 
order you have placed with NTlS, please Call our Customer Service 
Department at (703) 487-4660. 

About NTIS 

NTIS Collects scientific, technical, engineering, and business related 
information — then organizes, maintains, and disseminates that 
information in a variety of formats — from microfiche to online services. 
The NTIS collection of nearly 3 million titles includes reports describing 
research Conducted of sponsored by federal agencies and their 
contractors; statistical and business information; U.S. military 
publications; audiovisual products; computer software and electronic 
databases developed by federal agencies; training tools; and technical 
reports prepared by research organizations worldwide. Approximately 
100,000 new titles are added and indexed into the NTlS collection 
annually. 


For more information about NTIS products and services, call NTIS 
at (703) 487-4650 and request the free NTIS Catalog of Products 
and Services, PR-827LPG, or visit the NTIS Web site 
http://www.ntls.gov. 


NTIS 

Your indispensable resource for government-sponsored 
information— U.S. and worldwide 



